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Abstract 

Bipartite networks are composed of two types of nodes and there are no links between nodes of the same 
type. Thus the study of epidemic spread and control on such networks is relevant to sexually transmitted diseases 
(STDs). When entire populations of two types cannot be immunized and the effect of immunization is not per- 
fect, we have to consider the targeted immunization with immunization rates. We derive the epidemic thresholds 
of SIR and SIS models with immunization and illustrate the results with STDs on heterosexual contact networks. 
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1. Introduction 

Bipartite networks or graphs are composed of two 
types of nodes (or vertices) and there are no links (or 
edges) between nodes of the same type. Links are al- 
lowed to connect nodes of different types. Thus the 
study of epidemic spread and control on such networks is 
important for sexually transmitted diseases (STDs). A 
typical method to prevent or control an epidemic spread 
is immunization such as vaccination [2]. However, it is 
generally impossible to perfectly immunize entire popu- 
lations of two types, so objects or target sets of immu- 
nization and immunization rates must be introduced. 
Immunization rates are parameters representing various 
effects of immunity. 

We consider the SIR and SIS models with immuniza- 
tion targets and rates on bipartite networks. Based on 
[10, 11] two (partial) immunization policies are treated: 
"imperfect immunization to targeted nodes" and "insuf- 
ficient immunity" . 

The degree of a node is the number of links ema- 
nating from it. Using the probability distributions of 
degrees of both types, we derive the critical infection 
rates or the thresholds, above which a disease spreads 
in a bipartite network and below which it dies out. 

Due to [3, 4], sexual contact networks relevant to 
STDs follow power-law degree distributions; 

piku) oc kM~^^ and p{kY) oc fcp"'"', 

regardless of reported countries (Sweden, U.K., Burkina 
Faso, etc.) and periods (12 months and the life span). 
Here indices M and F mean two populations of males 
and females, respectively, and the exponents 7m and 7f 
are positive constants and typically satisfy 

2 < 7m < 3 < 7f. 

This paper deals with the SIR and SIS models with 
immunization, but the models involving the target sets 
and immunization rates are restricted to the SIS mod- 
els, because the derivation of the thresholds for the SIR 
models becomes a little lengthy. 

Finally, based on the thresholds, we apply immu- 
nization strategies to STDs on heterosexual contact 



networks as above. The immunization strategies should 
be implemented on nodes with high degrees in one pop- 
ulation and on those with low degrees in the other. 



2. The SIR model without immunization 

We distinguish two types of nodes of a bipartite net- 
work by indices 1 and 2. Hence all links connect nodes 
between types 1 and 2. In the SIR model, nodes of the 
network are divided into the following three groups re- 
garding infection states of a disease ([2], [6, Chap. 10]): 
Susceptible (S), Infected (I) and Removed (R). A sus- 
ceptible node (S-node) of type i = 1,2 becomes an in- 
fected node (I-node) at a rate Ai > 0. The parameters 
\i are the infection rates, for which we will derive the 
critical values for an epidemic outbreak. The disease 
can be passed from I-nodes to S-nodes of different types 
through links on bipartite networks. A recovered node 
(R-node) has either recovered from the disease or died 
and so they cannot pass the disease to others. An I-node 
becomes an R-node at a rate 5i {5i > 0) for type i. We 
can assume Si = 1, without loss of generality. 

Let p{ki) and p(k2) be the probability distributions 
of nodes with degree ki in population 1 and of nodes 
with degree ^2 in population 2, respectively. The aver- 
ages or moments are defined by 



for type i. 

Within population of each type i, the densities of 
S-, I-, R- nodes with degree ki at time t are denoted by 
variables S'fe . {t) , pfe . {t) , _Rfe . {t) , respectively. Since there 
are no rewiring of links, we have 

SkAt) + Pk^t) + Rk,{t) = 1. 

Following the dynamical mean-field approach ([5, 8, 
9]), we see that the spreading process in population 1 
on a bipartite network can be described by the system 



1 



of differential equations: 
dSk, 



dt 

dpki 
dt 

dRki 
dt 



= -XihSkAt)02{t), (1) 
= \ihSkAt)02{t) - 6ipkAt), (2) 
= SiPkAt)- (3) 



The term XikiSki{t)92{t) in (1) and (2) indicates the 
fraction of newly infected nodes through ki links, while 
02{t) is the probability of contact with I-nodes of pop- 
ulation 2 from which the disease spreads. Hence 02{t) 
can be written as 



^ _fe2 



(4) 



By the bipartite nature, analogous equations in pop- 
ulation 2 to Eqs. (1)- (4) arc valid under the interchange 
with 1 and 2. We indicate those equations by Eqs. (!')- 
(4'), respectively. The initial conditions of Ski set 
by 

5fe,(0) = l (i = l,2), (5) 

meaning that almost all nodes arc S-nodes at first. 

To begin with we solve Eqs. (1) and (!') under the 
initial conditions (5); 

5fe,(i) = e-^i'=i^^W, 5fe,(t) = e-^^'==^iW, (6) 

where 

Mt) = I 02{T)dT, Mt) = I 0l{T)dT. 

Jo Jo 

Using (3'), (4) and i?fe2(0) = for all k2, the auxiliary 
function (j)2{t) has an expression; 

1 /■* 
Mt) = Trryi / k2p{k2)pkA'r)dT 

We derive the differential equation for ^2{t)- Using 
(6), it follows that 



= E k2P{k2){l - Sk2 {t) - Rk2 it)) 



1 - 77^ V fep(fc2)e-^^'==^^W _ s^Mt)- 
{k2) ^ 



We are concerned with a steady state of the epidemic 
outbreak, in which one has the limits 



$i = lim ^i(t) [i = 1,2), 

t— >-oo 



together with the conditions 

iit) 



lim 

t^-cx) dt 



= (i = l,2). 



Using these relations in Eq. (7) and then interchanging 
1 and 2, we get the coupled equations for $j as follows: 



^2*2 = 1 - E k2P{k2)e-^''''^' , 



(8) 



= 1 - -L V fcip(fci)e-^i'=i*^ (9) 



Note that $2 = if $i = and vice versa, so the 

epidemic outbreak occurs when $1 > 0. The right hand 
sides of Eqs. (8) and (9) are monotone increasing and 
concave functions of $1 and $2, respectively. Therefore, 
the function 



which is derived from (9), is also a monotone increasing 
and concave function of <l>i via Eq. (8), and its value at 
$1 = 1 is less than 1 from the assumption di = 1. Thus 
the condition for the outbreak is given by 



d 1 



Remarking the relation 



d^i 



52{k2) ' 



it follows that the critical infection rates of Ai and A2, 
or thresholds, satisfy 



(10) 



This condition coincides with Eq. (5) in [3], which was 
obtained for the SIS model of STDs without immuniza- 
tion. Also in [7] a similar expression was derived using 
the generating function methodology. 

Since the SIS model does not contain R-nodes, it is 
easier to deal with the epidemic spreading with immu- 
nization within the framework introduced in Section 1. 



3. The SIS models with immunization 
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In this section we consider the spreading and control 
of a disease based on SIS models with target sets and 
iniiimnization rates. In these models R-nodes are absent 
and nodes that are recovered from the disease instantly 
become S-nodes again. 

A typical method to prevent or control an epidemic 
spreading is immunization such as vaccination [2] . Since 
it is very difficult or impossible to perfectly immunize an 
entire population, the object or target of immunization 
must be prescribed. Within population i {i = 1,2) we 
prescribe the object or target set, Tj, for immunization. 

We assume that the target sets are characterized in 
terms of degrees fcj. The notation fcj e Tj means that 
the population of all nodes with degree ki is an object of 
immunization. In case of Ti = {fci|A:i > i^i}, for exam- 
ple, the population of all nodes with degrees exceeding 
a size K\ is collectively an object of immunization. The 
notation Ti indicates the complement of Ti, meaning 
that Ti is not an object of immunization. The summa- 
tion restricted to all fc, in Tj will be denoted by X^^.. , 
and similarly for Ti by X^t ■ Furthermore, the averages 
of ki and kf over Tj are denoted by (fci)^. and (A;?)^. : 

k,eT, Ti 



Similarly the averages (fcj)^ and over T, are also 

defined. 

Two SIS models below involve immunization targets 
Tj and rates ccj (0 < aj < 1). The rates Ui are assumed 
to be constants, although they can be dependent on ki. 
The condition ai = 1 implies the perfect immunization 
for Ti, while aj = means no immunization. The vari- 
ables Sk^{t), pk.{t) indicate the densities of S-, I-nodes 
with degree ki at time t, as before. 



A. Imperfect immunization to targeted nodes 

This is the case where all of target nodes in Ti may 
not be immunized, because some arc overlooked or hid- 
den. The immunization rates ai (0 < < 1) represent 
the effectiveness of immunization such as vaccination 
coverage. 

Since immunized nodes are no longer S-nodes, we 
have 



1 - ai - Pk,{t)i 



for ki €Ti, 
for ki G Ti. 



Therefore, Eqs. (l)-(3) in Section 2 are replaced by the 
differential equations for population 1 



dt 



AiA:i(l — ai 



PkAt))02{t)-5ipkAt), 

if ki e Ti, 



Xikiii - pkAt)W2{t) - SiPkAt), _ 

if ki G Ti, 



(11) 



where 02 {t) is the same probability as (4) . Similar equa- 
tions for pk^it) follow from Eqs. (1') (3') and 9i{f) of 
(40- 

At the steady state, as in Section 2, we will have the 
conditions 

lin,%W=0, lim%i^=0 
t->oo dt t-foo dt 

for all ki and k2, and the limits 



91 = lim 0i{t), 62 = lim 6*2 (t). 



So we get from (11), 



hm pfe, {t) 

t— )-CSO 



f (l-ai)Aifcie2/(^i+Aifcie2), 
if fci eTi, 

Aifcie2/(^i + Aifcie2), 

if ki eTi. 



Substituting these into (4') in the limit, we have the 
equation for ©1 as follows: 



61 



1 - q;i)Ai 



Ti 



klp{ki)Q2 
5x + \ikiQ2 

klv{ki)Q2 
5i+\iki(c)2 



(12) 



By interchanging 1 and 2, we also have 



T2 



klp{k2)Qi 

^62 + A2fc26l 

^2 + A2A;2©i/' 



(13) 



If the coupled Eqs. (12), (13) have solutions 0i > 
and ©2 > 0, then an endemic outbreak occurs. The 
right hand sides of Eqs. (12) and (13) are monotone 
increasing and concave functions of ©2 and ©1, respec- 
tively. Hence the right hand side of (12) becomes a 
monotone increasing and concave function of ©1 via Eq. 
(13), and furthermore its value at ©1 = 1 is less than 
1. Therefore, we see that the thresholds of Ai and A2 
satisfy 



d 
d© 



+Eai 



(l-ai)Ai 

fc?p(fcl)©2 

61 + \1k1Q2 



klp{ki)Q2 
5i + AiA;i©2 



61=0 



1. 



Noting 



A2 



rf©2(0) _ 

rf©i ~ ^2(A;2) 



((i-«2)(fci>T. + (^i)fJ. 
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and using the identities 

we see that at the thresholds of Ai and A2, the relation 



AiAo = 



((fcf)-ai(fcf)^J((fc|)-a2(fc2%) 



(14) 



holds. In the case of no immunization {ai = 0:2 = 0) it 
coincides with (10). 

B. Insufficient immunity 

The parameters of this model represent the lev- 
els of immunity. The higher they are, the less infected 
S-nodes are. Thus we have 

c m_ / (l-ai)(l-Pfe,(t)), forfc, eT„ 
^kAt)-<^ l-pfc,(t), forfciGTi. 



So instead of (11) we obtain 



{ Aifci(l - ai)(l - PkAt))02{t) - Sipkdt), 

if ki e Ti, 



dt 



if fci e Ti. 



Similar equations for (^) ^-Iso hold by interchanging 
1 and 2. This model is the same as the one treated in 
[1] for unipartitc networks. Repeating the above calcu- 
lations in A, it follows that the condition for thresholds 
takes the same form as (14). 

As for the condition of the thresholds for SIR mod- 
els with immimization policies A and B, we are able to 
derive it by dividing each population i into two parts Tj 
and Ti as in [10], and using the procedures in Section 
2. Then we see that the same condition (14) is also 
obtained for on bipartite networks. Thus we conclude 
that the SIR and SIS models with immunization rates 
all have the same condition of thresholds as expressed 
in (14). 

Furthermore, populations 1 and 2 may adopt differ- 
ent immunization policies. Even if population 1 adopts 
the policy A, while population 2 adopts the policy B, 
we get the same condition of thresholds again. 

4. Immunization strategies 

Suppose that the degree distributions of a bipartite 
network follow power-laws; 

p{ki) (X ki^ and p{k2) oc k2~^^ 

for types 1 and 2, respectively. Moreover, the two expo- 
nents satisfy 

2 < 71 < 3 < 72, 



as in most sexual contact networks [3]. Other cases such 
as 

2<7ii72<3, 3<7i>72, 

can be discussed in similar ways. 

By (14) the geometric mean of the two critical 
thresholds 



Ac = 



Sld2{ki){k2) 



m) 



Ot2{kl)j.^) 



provides an overall threshold of the epidemic spreading. 
In order to control the spreading of STDs, we try to 
raise the value of Ac by employing two immunization 
strategies. It is desirable to implement both strategies 
simultaneously. 

Under the condition 2 < 71 < 3 the mean square 
(A;^) tends to infinity in the limit of infinite population, 
although real networks have finite sizes. First of all we 
must set Q!i = 1 for population 1. For, if ai < 1, then 

{kl) - at{kl)^^ > (1 - at){kl) 00 

in the limit of infinite population, no matter how Ti is 
chosen. So by setting ai = 1, we have 



(^i)ti 



<Kl 



for a target set Ti = {fci|fci > Ki\. Therefore, an im- 
munization strategy to population 1 is the perfect im- 
munization to nodes with an upper half of degrees as in 
[9]. 

On the contrary, (fc|) is finite by 72 > 3. This implies 
that nodes with high degrees make little contribution to 
and hence they are negligible. Thus one should take 
T2 — {^21^2 < K2}^ for some positive K2, bs, & target 
set of population 2. Then we have 



{kl) - cx2{kl)^^ = {kl)f^ + (1 - a2)(fci>, 



0, 



when a2 1 and K2 becomes larger. Hence the im- 
munization strategy to population 2 is concentrated on 
nodes with a lower half of degrees. Since is also 

finite, a2 need not be one, and as for K2 generally a 
small size is enough. 
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